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(57) Abstract: The invention concerns a method 
for processing an image acquired through a guide 
consisting of a plurality of optical fibers. The in- 
vention is characterized in that it consists, for each 
optical fiber, in isolating on the acquired image 
a zone corresponding to said optical fiber, in lo- 
cally processing each zone individually to correct 
the photon flux detected in each optical fiber, then 
in reconstructing the acquired image by eliminat- 
ing the pattern caused by the optical fiber. The 
method also comprises a sampling process for ob- 
taining, for each optical fiber and fix)m a sampling 
image, a sample injection rate which can be used 
for reconstructing the acquired images. The in- 
vention also includes a prior step which consists 
in detecting the fibers from a target. 

(57) Abrege : U invention conceme un procede 
de traitement d' image acquise au moyen d'un 
guide constitue par une pluralite de fibres 
optiques. Selon 1' invention, pour chaque fibre 
optique, on isole sur 1' image acquise une zone 
correspondante a cette fibre optique, on traite 
localement chaque zone de fa^on individuelle 
pour corriger le flux de photon d^tect^ dans 
chaque fibre optique, puis on reconstruit T image 
acquise en 61iminant le motif du aux fibres 
optiques. Le procede comprend egalement un 
processus d'etalonnage permettant 
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d'obtenir, pour chaque fibre optique et k partir d'une image d*etalonnage, un taux d' injection etalon qui va servir a la reconstruction 
des images acquises. Le proc6d6 comprend egalement en premier lieu une 6tape de detection des fibres k partir d*une image t^moin. 
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" Procede de trailsemen'b d'une image acquise au moyen d'un 
guide coxnpose d'une plurali1:e de fdLbres optl<iues." 



5 

La pr6sente invention se rapporte k un precede de traitement 
d'une image acquise au moyen d'un guide compose d'une pluralite 
de fibres optiques . Elle trouve une application particulidrement 
10 int^ressante dans le domaine de I'imagerie medicale. Toutefois 
1' invention est d'un cadre plus large puisqu'elle peut 
s'appliquer a tout domaine dans lequel on realise de I'imagerie 
au moyen d'un guide compose d'une pluralite de fibres optiques - 
Le guide d' image permet d'obtenir une image. Un tel dispositif 
15 permet de deporter le balayage laser et la source lumineuse et le 
recepteur loin de I'objet a observer. Par exemple dans un systeme 
a balayage laser avec obtention d'une image confocale, le guide 
d' image est un assemblage de plusieurs milliers de fibres 
optiques dont 1 ' arrangement spatial est identique en entree et en 
20 sortie. L ' observation d'un objet a travers ce guide pourrait etre 
compar§e a une observation via une grille , ^ cause de la perte 
d' information entre les fibres optiques. La visualisation est 
done g§nee a cause de la presence des fibres optiques : le motif 
des fibres optiques apparalt sur 1' image acquise. Cela impose un 
25 traitement sp6cifique de fagon k 61iminer ce motif et am^liorer 
la lisibilit6 de 1 'image - 

Dans les systemes actuels, ce traitement specif ique se 
limite a un filtrage lineaire de 1 ' image acquise. 

La presente invention a pour but de proposer un nouveau 
30 procede permettant de rendre les images acquises au moyen d'un 
guide multi-fibres optiques lisibles. 

Un autre but de 1' invention est de prendre en compte les 
effets parasites dus A I'appareil d' acquisition dans le 
traitement de 1' image acquise. 
35 On atteint au moins I'un des objectifs pr6cit6s avec un 

nouveau precede de traitement d' image acquise au moyen d'un guide 
constitu6 par une plurality de fibres optiques. Selon 
1' invention, pour chaque fibre optique, on isole sur 1' image 
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acquise une zone correspondante a cette fibre optique, on traite 
localement chaque zone de fagon individuelle, puis on reconstruit 
1' image acquise en eliminant le motif da aux fibres optiques. 

Avec le proc6d6 selon 1' invention^ le fait d'isoler les 
5 fibres optiques sur 1' image, revient k isoler sur 1« image, la 
zone correspondant 4 chaque fibre. Lorsque les pixels 
representant 1' injection majoritaire en terme de surface (zone 
d' influence) dans une fibre sont isoles, on peut alors effectuer 
des traitements locaux sur chaque fibre optique. La faible 
10 diaphonie du guide et les reglages de 1' injection permettent de 
garantir le fait que le contenu informatif de chaque fibre ne 
depend pas de la fibre voisine, mais uniquement de la coherence 
spatiale de I'objet observe. L'appareil realisant 1 'acquisition 
d' image est r6gle de fagon k avoir suffisarament de pixels par 
15 fibre : on peut ainsi estimer pr6cis6ment 1 ' information qui est 
detectee par cette fibre et qui est r6partie sur les pixels 
representant la fibre optique. 

Ainsi, au lieu de traiter 1' image de fagon globale en 
r6alisant un simple filtrage lineaire comme dans I'art ant^rieur, 
20 dans le proced6 selon 1' invention on isole chaque fibre optique 
sur 1' image acquise et on traite 1 ' information detectee par 
chaque fibre optique. Avantageusement, l'appareil realisant 
1' acquisition d' image est controle pour garantir les conditions 
minimales d'efficacite du precede selon 1' invention. Pour ce 
25 faire, on peut modifier le taux d' echantillonnage, la qualite 
d' injection dans les fibres optiques, et le reglage de la chaine 
de detection afin de garantir un profil de type "boite ^ oeufs", 
notamment sur 1 ' image temoin. 

Avantageusement, lorsque les fibres optiques sont isolees 
30 sur 1' image et lorsque chaque zone isol6e est traitee, de 
norabreuses applications peuvent etre envisagees telles que : 

la reconstruction d'une image sans le motif des fibres 
: les fibres g§nent la lisibilit6 et les traitements ulterieurs 
effeCtues sur 1' image; 
35 . le controle de la rugosity de la surface du guide, la 

rugosit6 va perturber le ph6nomSne d' injection pour le rendre 
spatialement variant; 

le recalage des images, ou la stabilisation de 
1' image; le motif des fibres emp§che tout recalage des images 
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entre elles^ la connaissance de la place des fibres et de 
1 ' information observ§e permettent de recaler les images; 

la super-r6solution : on pent utiliser des faibles 
mouvements dans 1 ' acquisition d'une sequence d' images pour re- 

5 echantillonner 1' image avec une periode spatiale plus faible, et 
done obtenir une meilleure resolution; 

la quantification des images : on peut extraire 
beaucoup plus facilement et pr^cisement des informations sur 
1' image sans le motif des fibres optiques, 

10 - le contrdle temporel des parametres internes de 

I'appareil d' acquisition : la connaissance de la place de chaque 
fibre et de leurs valeurs optimales d' injection permet de 
controler I'usure du guide, et les variations de certains 
parametres optoelectroniques . 

15 Dans le precede selon 1' invention, pour isoler chaque zone, 

on peut appliquer un masque, correspondant au motif des fibres, 
sur 1' image acquise. Ce masque, correspondant a une image des 
composantes connexes representant chaque fibre ^ est obtenu au 
cours d'une etape de detection des fibres ^ partir d'une image 

20 temoin . 

L'iiaage temoin est une image permettant de bien distinguer 
les fibres optiques les unes des autres. Elle peut etre issue de 
1 ' observation d'un miroir, d^un milieu diffusant homog^ne, d'un 
milieu fluorescent homogene, elle peut aussi provenir de la 
25 retro-diffusion propre a 1' inter ieur du paquet de fibres 
optiques. Mais elle peut encore etre 1' image acquise. A la sortie 
de la detection, on obtient done une image des composantes 
connexes (des segments) representant chaque fibre optique. Chaque 
niveau de gris represente un index unique designant une fibre 
30 optique dans le guide. 

Selon 1' invention, 1' etape de detection des fibres peut 
comprendre les etapes suivantes : 

pre-filtrage de 1' image temoin, 

segmentation par region, en utilisant I'algorithme de 
35 "ligne de partage des eaux" LPE, 

correction de segments pr^sentant une surface 
anormalement grande, et 
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- correction de segments pr6sentant une surface 

anormalement petite. 
Les deux etapes de corrections sont interchangeables, et 
elles peuvent Stre r6alisees de fagon iterative. 
5 Avantageusement,. I'^tape de pre-filtrage peut comprendre une 

etape d'ouverture morphologique suivie d'une 6tape d' inversion 
d* image. Avec I'ouverture morphologique numerique, on cherche a 
eliminer les maxima parasites situes sur les fibres optiques, 
C'est un pr§-traitement classique de I'algorithme de "ligne de 
10 partage des eaux" LPE effectue lors de la segmentation par 
region. 

A la sortie du pre-filtrage, on obtient une image des fibres 
optiques filtrees de leurs maxima locaux, et lissees au niveau 
des zones inter-fibres, 
15 L' etape d' inversion d' image peut etre precedee d'une 6tape 

de diffusion anisotrope de type scalaire. 

Selon 1' invention, le pr6-filtrage peut comprendre en outre 
une etape au cours de laquelle on realise une interpolation au 
plus proche voisin pour doubler la taille de 1' image en vertical 
20 et en horizontal. 

Avec cette §tape d' interpolation, on cherche ^' simuler des 
Elements structurant de morphologie mathematique avec un rayon 
inf^rieur a un. L' image est doublee pour que I'ouverture 
morphologique qui suit ne touche pas aux maximum isoles, mais 
25 seulement ceux qui sont 8-connexe, mais non 4-connexe (voisins 
par une diagonale) . L'interet est de faire une selection des 
maxima elimines par I'ouverture. 

Par ailleurs, en presence d'une pluralite d' images 
d' acquisition, le pre-filtrage peut comprendre en outre une etape 
30 de f iltrage temporel . 

Avantageusement, le traitement local de chaque zone peut 
consister ^ calculer le flux de photon detecte pour chaque zone 
(correspondant ^ une fibre donnee) de 1' image acquise, et a 
corriger le biais sur chaque valeur de flux ainsi calculee. 
35 De pr6f6rence, on realise le calcul de flux au moyen d'un 

estimateur du maximum de vraiserablance calcul6 sur un profil 
sp§cifique d' injection de chaque fibre. Plus precisement, on peut 
utiliser 1' estimateur du maximum de vraisemblance sur la 
distribution d'amplitude du profil specif ique d' injection dans 
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chaque fibre optique. Le profil est une courbe repr6sentant le 
taux d' injection en fonction de l'§loigneiaent de la lumiere par 
rapport au centre de la coupe transversale k l'extr6niit§ de la 
fibre optique. Souvent ce profil est modelise par une gaussienne. 
5 Suivant un mode de mise en oeuvre de 1' invention, en 

appliquant egalement le masque sur une image representant un fond 
parasite, on calcule egalement le flux de photon detecte pour 
chaque zone de 1' image de fond, et on soustrait a chaque valeur 
de flux de chaque zone de 1' image acquise, la valeur de flux de 
10 chaque zone de 1' image de fond correspondante, et on realise la 
correction de biais sur le resultat de cette soustraction. 

L' image de fond peut etre les reflexions parasites sur les 
systemes optiques de I'appareil d' acquisition, et done y compris 
sur la sortie du guide d' image, mais elle peut egalement etre le 
15 decalage, le bruit electronique, de la chaine de num§risation de 
I'appareil d' acquisition. Le d6calage correspond au terme 
communement appel§ "offset". Si 1' offset est dominant sur 
1' image, on ne peut pas obtenir le fond simplement en retirant 
1' image, car 1* offset depend du contenu, et n'est done plus le 
20 meme. Dans ce cas, on utilise un quantile de 1 ' histogramme pour 
I'es timer. L' histogramme est celui de 1' image acquise lors de la 
mesure en temps reel, et celui d'une image d'etalonnage lors 
d'une etape d'6talonnage comme on le verra ci-dessous. 

Le fond parasite peut provenir du fond de 1' image ou d'un 
25 offset de la chaine de detection. 

Selon 1? invention, la correction de biais peut consister a 
separer spatialement les fibres en differents blocs, a estimer la 
valeur de biais dans chaque bloc, a interpoler les valeurs des 
biais de fagon ^ obtenir une valeur de biais pour chaque fibre, 
30 et a diviser, pour chaque zone, la valeur de flux obtenue A 
1' etape precedente par la valeur de biais correspondante ainsi 
obtenue . 

La reconstruction de 1' image acquise peut faire intervenir 
une etape de calibrage pour calibrer le flux de 1' image acquise, 
35 apr^s traitement local, et une 6tape de reconstruction mosaique. 
On peut utiliser d'autres types de reconstruction, telles que par 
interpolation ou avec des bases de fonctions radiales. 
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Par ailleurs, la presente invention peut etre mise en oeuvre 
sans les 6tapes d'6talonnage et de calibrage. Dans ce cas, 
1' image temoin peut etre 1' image acquise. 

Pour le calibrage et pour chaque zone de 1' image acquise, on 
5 peut diviser la valeur du flux obtenue aprSs traitement local par 
une valeur de flux obtenue A la suite d'une etape d' 6talonnage . 
Cette operation de division permet de compenser les mauvaises 
injections dans certaines fibres optiques. 

Selon une caracteristique avantageuse de 1' invention, 
10 1* etape d'etalonnage consiste a : 

isoler chaque zone d'une image d'etalonnage en appliquant 
le masque, correspondant au motif des fibres, sur cette 
image d'etalonnage, 
- calculer le flux de photon detecte pour chaque zone de 
15 1' image d'etalonnage, et 

corriger le biais sur chaque valeur de flux ainsi 
calculee. 

L' image obtenue a 1' issue de I'etalonnage pourra servir 
d'etalon a 1' image acquise de fagon a obtenir une image acquise 
20 pour laquelle toutes les fibres optiques du guide auraient §te 
inject§es de la m§me mani^re, 

De preference, on realise le calcul de flux au moyen d'un 
estimateur du maximum de vraisemblance calcule sur le profil 
specif ique d' injection de chaque fibre. Par ailleurs, en 
25 appliquant egalement le masque sur une image repr^sentant un fond 
parasite, on peut calculer le flux de photon detects pour chaque 
zone de 1' image de fond, on peut soustraire a chaque valeur de 
flux de chaque zone de 1' image d'etalonnage, la valeur de flux de 
chaque zone de 1 ' image de fond correspondante, et on peut 
30 realiser la correction de biais sur le resultat de cette 
sous tract ion. 

En d'autres termes, au cours de 1 ' etape d'etalonnage, on 
realise les memes operations que lors de la mesure, c'est-A-dire 
lors du traitement d'une image acquise en temps r6el. Seulement, 
35 d I'etalonnage on se sert d'une image d'etalonnage faisant 
largement apparaitre le motif des fibres optiques. L'etalonnage 
permet, apres detection des fibres sur 1' image temoin, de generer 
une image dans laquelle le taux d' injection va servir d'etalon 
lors de la mesure en temps reel. Au cours de la mesure, on 
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corrige 6galement le taux injection sur 1' image acquise, et on 
calibre le flux observ6 en fonction de 1' image etalon de fagon 4 
reconstruire une image acquise sans motif des fibres optiques. 

Par biais, on entend une composante basse frequence, cette 
5 composante pouvant provenir de di verses causes- 

La correction de biais peut etre effectu6e lors de 
I'etalonnage et lors de la mesure en temps reel. Pour le premier 
cas^ cela peut venir du fait que I'etalonnage se fait par exemple 
sur un miroir plan, et que la courbure de champ va reduire la 
10 qualite d' injection au retour sur les bords (qui sont 
defocalises) . Sur I'objet de mesure, ou dans un milieu diffusant 
homogene, 1' injection reste moins bonne sur les bords, et cela se 
traduit par un biais trcs - similaire au premier quant a sa forme. 
Le biais peut egalement provenir d'un probl^me de vignettage. 
15 Generalement, le biais a une symetrie quasi circulaire. 
L' estimation du biais se fait en divisant 1' image en NxN blocs de 
taille fixe, puis en estimant le biais sur chaque bloc. Pour cela 
il faut considerer la nature de I'objet observe. Dans le cas d'un 
objet homogfene, le biais peut §tre acquis en prenant la valeur 
20 moyenne ou mediane sur le bloc. Quand il y a un objet, .il faut 
savoir si cet objet est plus sombre ou plus clair que le reste de 
1' image. A titre d' exemple, le biais peut %tre multiplicatif , et 
on prend done plutot un operateur de moyenne ou de mediane (par 
rapport & un max ou min pour un biais additif ) - On obtient alors 
25 une image de taille NxN qui est utilisee, apres interpolation, 
pour trouver la valeur du biais vue par chaque fibre. On peut 
utiliser une interpolation bilineaire par exemple. 

La reconstruction mosaique peut consister a repartir sur 
toute la surface de chaque zone de 1' image acquise, la valeur de 
30 flux de chaque zone obtenue a la suite de l'6tape de calibrage. 
On peut ensuite realiser un filtrage recursif passe-bas de fagon 
a lisser 1' image acquise reconstruite. 

Selon un mode de mise en oeuvre avantageux de 1' invention, 
1' image t6moin et 1' image d'etalonnage sont identiques . 
35 D'autres avantages et caracteristiques de 1' invention 

apparaitront a I'examen de la description detaillee d'un mode de 
mise en oeuvre nullement limitatif, et des des sins annex6s, sur 
lesquels : 
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La figure 1 est une vue sch6iaatique globale du proc6de de 
traitement d'' image selon 1' invention ; 

La figure 2 est un organigramrae dStaillant les 
principales etapes d'un processus d' 6talonnage selon 
5 1' invention ; 

La figure 3 est un organigramme detaillant les 
principales etapes d'un processus de mesure selon 
1' invention ; et 

La figure 4 est un organigramme detaillant les etapes 
10 finales prenant en compte les processus d' 6talonnage et 

de mesure pour la reconstruction d'une image acquise sans 
fibres optiques apparentes selon 1' invention. 
Sur la figure 1 est represents un schema global du precede 
selon 1' invention. 
15 Le procede peut etre scinde en quatre parties : une partie 1 

concernant la detection des fibres, une partie 2 d'' etalonnage, 
une partie 3 de mesure, et une partie 4 de reconstruction. Les 
parties 1 et 2 correspondent au diagramme de la figure 2, alors 
que la partie 3 correspond au diagramme de la figure 3, et la 
20 partie 4 correspond au diagramme de la figure 4. 

Sur la figure 1, selon 1' invention, lorsqu'on desire 
acquerir une s§rie d' image, on realise d'abord une etape 
d' Stalonnage . Pour ce faire, on considere une image temoin 5 a 
laquelle on fait subir une §tape de detection des fibres 6 de 
25 faq:on a obtenir une image 7 des composantes connexes reprSsentant 
chaque fibre. L' image temoin 5 est une image acquise au moyen 
d'un systeme comprenant un guide constitue d'une pluralite de 
fibres optiques de 10000 a 30000 par exemple. L' image temoin 5 
est obtenue de telle sorte que I'on distingue le motif des fibres 
30 optiques, c'est a dire de type "boite a oeuf" : sur le profil de 
1' image, une fibre se traduisant par une petite montagne entouree 
de col et de vallee. L' image temoin 5 subit une operation de 
detection des fibres de fagon 3i obtenir une sorte de masque 
representant le motif des fibres optiques. Ce masque est 1' image 
35 7 des composantes connexes repr6sentant chaque fibre. A chaque 
niveau de gris represente un index unique designant une fibre 
optique dans le guide. 
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L' image 7 est ensuite utilis§e pour l'6talonnage 2 du 
systeme d' acquisition d' image. L'6talonnage a pour but de 
determiner une image des taux d' injection de photons fibre k 
fibre. Cette 6tape d'etalonnage est n6cessaire dans la mesure ou 
5 chaque fibre pr6sente des propri6tes physiques 16g^rement 
diff6rentes des autres fibres. II existe done une certaine 
disparity concernant la capacity de chaque fibre optique k 
v6hiculer un mSme flux de photon. 

Pour r^aliser I'etalonnage 2, on considere une image 8, ci- 
10 dessous appel6e image miroir, obtenue en plagant le miroir devant 
le syst6me optique du guide d' image. Cette image peut egalement 
etre celle d'un milieu diffusant homog^ne, d'un milieu 
fluorescent homogene, ou de la r6tro-dif fusion propre k 
I'int^rieur du paquet de fibres optique. Cette image 8 peut 
15 egalement etre la m§me image utilisee en 5, c' est A dire 1' image 
t^moin. On utilise le masque 7 pour determiner le flux de photon 
de 1' image miroir 8 vue par chaque fibre optique au cours de 
I'etape 9. De fagon facultative, on peut egalement, determiner le 
flux de photon vu par chaque fibre au cours de l'6tape 11 pour 
20 une image de fond parasite 10. Cette image 10 peut correspondre k 
des reflexions parasites • sur les optiques du systeme 
d' acquisition, mais aussi a 1' offset et/ou le bruit eiectronique 
da h la chaine de numerisation du systems d' acquisition. Le 
calcul de flux de I'etape 11 fait egalement intervenir le masque 
25 7 de fagon k identifier la zone correspondante k chaque fibre. A 
I'etape 12, pour chaque fibre optique, on soustrait la valeur de 
photon de 1' image miroir 8 par la valeur de flux de photon de 
1' image de fond 10. A I'etape 13 on estime que pour chaque fibre 
optique, la difference obtenue k I'etape 12 correspond au taux 
30 d' injection etalon pour chaque fibre, (etape 13). 

A I'etape 14, on realise une correction de biais sur 1' image 
de I'etape 13. L' image issue de I'etape 14 est done une image 
presentant pour chaque zone correspondant k une fibre optique, 
une valeur de flux de photon etalon et corrige. Cette image issue 
35 de I'etape 14 va servir de reference pour une serie d' images 
acquises en temps reel par le systeme d' acquisition. Les images 
acquises et traitees en temps reel subissent le traitement 
illustre dans les parties 3 et 4 . 
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La partie 3 de mesure regoit une image acquise 15/ 
typiquement 1' image d'un objet de mesure. Comme realise dans la 
partie 6talonnage 2^ on calcule 6galement ici A I'etape 18 le 
flux de photon vu par chaque fibre en ce qui concerne 1' image 
5 acquise 15. Pour ce faire, on utilise le masque 7 de fagon a 
identifier sur 1' image acquise 15, la zone correspondant k chaque 
fibre optique. De fagon facultative, on considere une image de 
fond parasite 16 comme precedeiranent, qui peut etre une image, 
reelle, c' est ^ dire correspondant au fond de 1' image acquise 15 
10 ou alors une image estim^e correspondant au parasite du systeme 
d'^ acquisition. Cette image de fond 16 subit egalement une etape 
17 au moyen du masque 7 de fagon a determiner le flux de photon 
vu par chaque fibre. A I'etape 19 on realise une soustraction. Le 
resultat de la soustraction est une image presentant pour chaque 
15 zone correspondant a une fibre optique donn6e, un flux de photon 
utile- L'6tape 21 est une etape facultative au cours de laquelle 
on realise une correction de biais sur 1' image 20. 

La partie 4 de reconstruction regoit d'une part 1' image 20 
corrigee et d'une autre part 1' image 13 corrig6e (de-biaisee) de 
20 fagon A effectuer une operation de calibrage 22 en divisant les 
flux de 1' objet observ6 {616ment issu de I'etape 21) par les flux 
etalons (616ment issu de I'etape 14). A I'etape 22 on realise 
6galement une reconstruction de fagon k obtenir une image 
reconstruite 23 sans motif des fibres optiques apparent. 
25 Sur la figure 2 on voit plus en detail le processus 1 de 

detection des fibres et le processus 2 d' etalonnage . L' operation 
6 de detection des fibres fait intervenir quatre operations : 
un pre-f iltrage, 

une "ligne de partage des eaux" LPE correspondant a 
30 une segmentation par region ; 

une correction des segments ayant une surface 
anormalement grande ; et 

une correction des segments ayant une surface 
anormalement petite. 
35 Les deux operations de correction sont interchangeables, et 

on peut egalement les effectuer en boucle. 

L' operation de pre-f iltrage 61 regoit en entree 1' image 
temoin 5 et genfere une image des fibres filtrees de ces maxima 
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locaux^ et lissees au niveau des zones inter-fibres. Le pre- 
filtrage fait intervenir une operation d'ouverture morphologique , 
suivi event uellement d'une diffusion anisotrope de type scalaire, 
puis d'une inversion de 1' image. Lorsqu'on est en presence de 
5 plusieurs images d'un meme objet fixe, on peut effectuer un 
filtrage temporel des images. On peut egalement effectuer une 
interpolation au plus proche voisin pour doubler la taille de 
1' image en vertical et horizontal. 

L' image generee par le pre-f iltrage 61 subit ensuite une 
10 operation de ligne de partage des eaux 62, permettant I'obtention 
d'une image des composantes connexes des fibres detectees. 
L' operation de ligne de partage des eaux, de type conventionnel, 
va permettre de reperer des segments se trouvant au bord de 
1' image de fagon a les retirer du r6sultat final. On calcule 
15 ensuite les caract6ristiques de ces segments obtenus (taille 
moyenne, ecart type, voisinage) . L' image issue de 1' operation 62 
va subir deux corrections 63 et 64 successives. L' operation 63 
est une correction des segments ayant une surface anormalement 
grande. Pour ce faire, on selectionne les segments qui sont a la 
20 fois trop gros par rapport a une taille moyenne et qui ont trop 
de voisins par rapport a une surface normale des fibres. Ces 
segments sont re-segmentes avec une ligne de partage des eaux 
soit sur 1' image d'origin^, soit sur une image de carte de 
distance a I'interieur des segments detectes (dans un segment, 
25 distance entre chaque pixel et le bord du segment) . Puis, on 
calcule les caracteristiques des segments obtenus (taille 
moyenne, ecart type des tallies) . 

En 64, on corrige les segments ayant une surface 
anormalement petite. On calcule aussi le graphe d'adjacence des 
30 segments, puis on decide quelles fibres doivent etre 
obligatoirement fusionnees et celles qui sont juste candidates. 
Dans tous les cas, 1' ensemble des fusions possibles est 
1' ensemble des fusions avec chaque voisine. Pour les fibres 
obligatoirement fusionnees, on prend la fusion possible qui donne 
35 la valeur de compacit6 la plus petite. Pour les autres qui sont 
candidates seulement, on utilise successivement trois filtres 
pour eliminer les fusions qui donnent de mauvais resultats . Le 
premier f litre verifie que la taille apres fusion n' est pas trop 
grande. Le second f litre verifie que la compacit^ apres fusion ne 
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depasse pas une valeur maximale. Le dernier filtre verifie que la 
fusion ameliore la coiapacit§. S'il reste plusieurs fusions, on 
garde celles qui donnent un resultat de meilleure compacit§ (la 
plus petite) . Une fois toutes les fusions effectufees sur le 
5 graphe d'adjacence, on r6percute les r6sultats sur 1' image des 
composantes connexes, en sortie. Puis on calcule les 
caracteristiques des segments obtenus (taille moyenne, ecart type 
des tallies). La seconde correction 64 permet de generer 1' image 
des composantes connexes qui va servir de masque pour 
10 I'etalonnage 2 et la mesure 3. 

La partie droite de la figure 2 concerne un processus 
d'etalonnage 2 telle que represente sur la figure 1 mais de fagon 
simplifiee- En effet, sur la figure 2, les operations 
optionnelles 10 et 11 n' apparaissent pas. 
15 L'etape 91 de masquage, consiste a rep^rer sur 1' image 

miroir 8, la zone ou surface correspondant a chaque fibre optique 
du guide. L' image 7 des composantes connexes sert de masque. A 
I'etape 92, pour chaque fibre optique, on calcule le flux 
provenant de I'objet observfe. Le flux est calcul6 grace k 
20 I'estimateur du maximum de vraisemblance d6termin6 sur le profil 
specif ique d'injection de chaque fibre optique. A I'etape 14, on 
realise une correction de biais comme on le verra plus en details 
sur la figure 3. A la sortie de I'etape 14, on obtient en 24 pour 
chaque fibre optique, une valeur 6talon du taux d' injection de 
25 photons . 

Sur la figure 3, on distingue le processus de mesure selon 
1' invention. Ce processus se fait en temps reel. Sur 1' image 
acquise 15, on repere la zone correspondante k chaque fibre 
optique en effectuant une operation de masquage 181 au moyen du 

30 masque 7. On calcule ensuite en 182, pour chaque fibre, le flux 
provenant de I'objet observe. Comme precedemment, le calcul est 
r6alis6 grace a I'estimateur du maximum de vraisemblance 
determine sur le profil sp6cifique d' injection de chaque fibre 
optique. On effectue les memes operations de masquage 171 et de 

35 calcul de flux 172 sur une image de fond parasite 16. 

A I'etape 19, pour chaque zone correspondant a une fibre 
optique, on soustrait le flux de fond (172) du flux de 1' image 
acquise (182) . 
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Si 1' image de fond n'est pas fournie et que ce fond existe, 
il faut le soustraire, on calcule alors un offset et/ou bruit 
61ectronique du systeme de mesure grace k un quantile sur 
1 ' histogramme de 1' image acquise 15. 

On peut ensuite r^aliser A l'6tape 21 une correction de 
biais sur 1' image issue de la soustraction 19. Dans ce cas, on 
separe spatialement en 211 les zones correspondant aux fibres en 
diff6rents blocs. Dans chacun de ces blocs, on calcul en 212 la 
valeur du biais grSce k un operateur donne. Puis en 213, on 
interpole les valeur s du biais afin d'obtenir une valeur pour 
chaque fibre. On realise alors en 214 une division de la valeur 
du flux vu par chaque fibre par la valeur du biais obtenu. 

A 1' issue de 1 ' etape 21, on gendre une image 25 reprfesentant 
le flux observe pour chaque fibre. 

L'ultime etape de reconstruction prend en compte 1' image du 
taux d' injection 6talon 24 et 1' image de fl\ax observe 25. 
L' operation de calibrage a pour but de compenser les pertes 
d' injection en egalisant le taux d' injection de toutes les fibres 
optiques de fa?on k avoir une image dont toutes les fibres ont 
6t§ injecte de la m§me fa9on. Pour cela, on divise en 221 1' image 
de flux observe 25 par 1' image du taux d' injection 24- 

On realise ensuite une reconstruction mosaique en 
rfepartissant k 1' etape 222 sur toute la surface correspondant k 
chaque fibre la valeur obteriue apres calibrage (division) . Afin 
de donner un aspect plus r^gulier, on peut effectuer un filtrage 
223 passe-bas gaussien par exemple. 

Bien sttr, 1' invention n'est pas limitee aux exemples qui 
viennent d'etre decrits et de nombreux am§nagements peuvent §tre 
apport6s a ces exemples sans sortir du cadre de 1' invention. 
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REVENDICATIONS 



10 



15 



20 



25 



30 



35 



1. Proced§ de traitement d' image acquise au moyen d'un guide 
constitu§ par une pluralite de fibres optiques, caract6rise en ce 
que, pour chaque fibre optique, on isole sur 1' image acquise une 
zone correspondante k cette fibre optique, on traite localement 
chaque zone de fagon individuelle, puis on reconstruit 1' image 
acquise en 61iminant le motif da aux fibres optiques. 

2. Proced^ salon la revendi cation 1, caracterise en ce que 
pour isoler chaque zone, on applique un masque, correspondant au 
motif des fibres, sur 1' image acquise. 

3. Precede selon la revendication 2, caract6ris6 en ce que 
le masque, correspondant A une image des composantes connexes 
repr6sentant chaque fibre, est obtenu au cours d'une 6tape de 
detection des fibres k partir d'une image temoin. 

4. Proc6d§ selon la revendication 3, caract6ris6 en ce que 
l'6tape de detection des fibres comprend les 6tapes suivantes : 

- pr6-filtrage de 1' image t6moin, 

- segmentation par region, 

- correction de segments presentant une surface 
anormalement grande, et 

- correction de segments presentant une surface 
anojnaalement petite . 

5. Proc6d§ selon la revendication 4, caracterise en ce que 
les deux etapes de corrections sont interchangeables . 

6. Precede selon la revendication 4 ou 5, caract6ris6 en ce 
que les deux stapes de corrections sont r6alis6es de fagon 
iterative . 

7. Proc6d6 selon I'une quelconque des revendications 4 k 6, 
caract§ris6 en ce que l'6tape de pr§-filtrage comprend une 6tape 
d'ouverture morphologique suivie d'une 6tape d' inversion d' image. 
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8. Proc6d6 selon la revendi cation 1, caracterise en ce que 
l'6tape d'inversion d* image est precedee d'une §tape de diffusion 
anisotrope de type scalaire. 

5 9. Proc6de selon I'une quelconque des revendi cations 4 a 8, 

caract6rise en ce que le pre-filtrage comprend en outre una etape 
au cours de laquelle on realise une interpolation au plus proche 
voisin pour doubler la taille de 1' image en vertical et en 
horizontal . 

10 

10. Precede selon I'une quelconque des revendications 4 a 9, 
caracterise en ce que, en presence d'une plurality d* images 
d' acquisition, le pre-filtrage comprend en outre une 6tape de 
filtrage temporel. 

15 

11. Procede selon I'une quelconque des revendications 
pr6c6dentes, caracteris6 en ce que le traitement local de chaque 
zone consiste ^ calculer le flux de photon detects pour chaque 
zone de 1' image acquise, et a corriger le biais sur chaque valeur 

20 de flux ainsi calculee. 

12. Proc6d6 selon la revendication 11, caracterise en ce 
qu'on realise le calcul de flux au moyen d'un estimateur du 
maximiam de vraisemblance calculi sur un profil specifique 

25 d' injection de chaque fibre. 

13. ProcedS selon la revendication 7, caracteris6 en ce 
qu'en appliquant egalement le masque sur une image representant 
un fond parasite, on calcule egalement le flux de photon detects 

30 pour chaque zone de 1 ' image de fond, et on soustrait ^ chaque 
valeur de flux de chaque zone de 1' image acquise, la valeur de 
flux de chaque zone de 1' image de fond correspondante, et on 
realise la correction de biais sur le resultat de cette 
soustraction. 

35 

14. Proced6 selon la revendication 13, caracterise en ce que 
le fond parasite provient du fond de 1' image. 
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15. Proc6de selon la revendication 13, caract§ris6 en ce que 
le fond parasite provient du calcul d'un decalage de la chaine de 
detection 

16. Proc6d6 selon I'une quelconque des revendi cat ions 11 a 
15, caract6rise en ce que la correction de biais consists a 
separer spatialement les fibres en differents blocs, a estimer la 
valeur de biais dans chaque bloc, a interpoler les valeurs des 
biais de fagon a obtenir une valeur de biais pour chaque fibre, 
et a diviser, pour chaque zone, la valeur de flux obtenue a 
I'etape precedente par la valeur de biais correspondante ainsi 
obtenue . 

17- Proc6d6 selon I'une quelconque des revendications 
15 pr6c6dentes, caracterise en ce que la reconstruction de 1* image 
acquise fait intervenir une etape de calibrage pour calibrer le 
flux de 1' image acquise, apres traitement local, et une 6tape de 
reconstruction mosaxque . 

20 18. Proced§ selon la revendication 17, caractSrise en ce que 

pour le calibrage et pour chaque zone de 1' image acquise, on 
divise la valeur du flux obtenue apres traitement local par une 
valeur de flux obtenue a la suite d'une etape d'etalonnage. 

19. Precede selon la revendication 18, caracterise en ce que 
1' etape d'etalonnage consiste a : 

isoler chaque zone d'une image d'etalonnage en appliquant le 
masque, correspondant au motif des fibres, sur cette- image 
d'etalonnage, 

calculer le flux de photon detecte pour chaque zone de 
1' image d'6talonnage, et 

corriger le biais sur chaque valeur de flux ainsi calculee. 

20. Precede selon la revenjdication 19, caracteris6 en ce 
35 qu'on realise le calcul de flux au mo yen d'un estimateur du 

maximum de vraisemblance calcule sur le profil sp6cifique 
d' injection de chaque fibre. 
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21. Proc6d6 selon la revendication 18 ou 19, caract6ris§ en 
ce qu'en appliquant 6galement le masque sur une image 
reprfesentant un fond parasite, on calcule egalement le flux de 
photon detecte pour chaque zone de 1* image de fond, on soustrait 
5 a chaque valeur de flux de chaque zone de 1' image d' etalonnage, 
la valeur de flux de chaque zone de 1' image de fond 
correspondante, et on realise la correction de biais sur le 
resultat de cette sous traction. 

10 22, Procede selon la revendication 21, caracterise en ce que 

le fond parasite provient du fond de 1' image. 

23 Procede selon la revendication 21, caracterise en ce que 
le fond parasite provient du calcul d'un d§calage et du bruit de 
15 la chalne de detection. 

24. Precede selon I'une quelconque des revendications 17 k 
23, caracterise en ce que la reconstruction mosalque consiste a 
r6partir sur toute la surface de chaque zone de 1' image acquise, 
20 la valeur de flux de chaque zone obtenue a la suite de I'etape de 
calibrage . 

25 Proced6 selon la revendication 24, caracterise en ce 
qu'on realise un filtrage passe-bas de fagon a rendre plus 
25 reguli^re 1' image acquise reconstruite . 

26. Procede selon I'une quelconque des revendications 3 ^ 
25, caracterise en ce que 1' image t6moin est une image obtenue en 
plagant un miroir face au guide . 

30 

27. Proc6d6 selon I'une quelconque des revendications 3 a 
25, caracterise en ce que 1' image temoin est une image obtenue a 
partir d'un milieu diffusant homogSne. 

35 28. Proc6de selon I'une quelconque des revendications 3 a 

25, caract6ris§ en ce que 1' image temoin est une image obtenue a 
partir d'un milieu fluorescent homogene 
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29. Proc6d6 selon l*une quelconque des revendications 3 A 
25r caract6ris6 en ce que 1' image t6iaoin est une image obtenue ^ 
partir de la retro-diffusion A l'int6rieur du paquet de fibres 
optiques constituant le guide. 

5 

30. Precede selon I'une quelconque des revendications 3 a 
25/ caracterise en ce que 1' image temoin est 1' image acquise. 

31. Precede selon I'une quelconque des revendications 19 a 
10 30 f caracterise en ce que 1 ' image temoin et 1' image d'etalonnage 

sont identiques- 

32. Appareil d' acquisition d' image au moyen d'lon guide 
constitue par une plurality de fibres optiques, et mettant en 

15 ceuvre un proc6d6 selon I'une quelconque des revendications 
prec§dentes, caracterise en ce que, pour chaque fibre optique, il 
comprend : 

des moyens pour isoler sur 1' image acquise une zone 
correspondante k cette fibre optique, 
20 - des moyens pour traiter localement chaque zone de fagon 

individuelle, et 

des moyens pour reconstruire 1 ' image acquise en eliminant le 
motif du aux fibres optiques. 

25 33. Appareil selon la revendication 32, caracterise en ce 

qu'il comprend des moyens pour modifier le taux 
d'^chantillonnage, la quality d' injection dans les fibres 
optiques, et le r^glage d'une chaine de detection afin de 
garantir un profil "boite a oeufs". 

30 

34. Application du proc6de de traitement d' image selon I'une 
quelconque des revendications pr6c6dentes pour I'un des domaines 
suivants : 

le controle de la rugosite de la surface du guide; 
35 - le recalage des images, ou la stabilisation de 1' images- 

la super-resolution d'une image acquise; 
la quantification des images; et 

le controle temporel des parametres internes de 1' appareil 
d ' acquisition . 
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